Quark spectral density and a strongly-coupled QGP 
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The maximum entropy method is used to compute the dressed-quark spectral density from the 
self- consistent numerical solution of a rainbow truncation of QCD's gap equation at temperatures 
above that for which chiral symmetry is restored. In addition to the normal and plasmino modes, the 
spectral function also exhibits an essentially nonperturbative zero mode for temperatures extending 
to 1.4 — 1.8-times the critical temperature, T c . In the neighbourhood of T c , this long- wavelength 
mode contains the bulk of the spectral strength and so long as this mode persists, the system may 
fairly be described as a strongly-coupled state of matter. 

PACS numbers: ll.10.Wx, 12.38.Mh, ll.15.Tk, 24.85.+p 



It is believed that a primordial state of matter has been 
recreated by the relativistic heavy- ion collider (RHIC) 
[1]. This substance appears to behave as a nearly-perfect 
fluid on some domain of temperature, T, above that re- 
quired for its creation, T c [2]. An ideal fluid has zero 
shear- viscosity: r\ = 0, and hence no resistance to the 
appearance and growth of transverse velocity gradients. 
A perfect fluid with near-zero viscosity is the best achiev- 
able approximation to that ideal. Graphene might pro- 
vide a room temperature example [3]. From Newton's 
law for viscous fluid flow; viz., (v/d) = (l/rj)(F/A), it 
is apparent that in near-perfect fluids a macroscopic ve- 
locity gradient is achieved from a microscopically small 
pressure. Strong interactions between particles consti- 
tuting the fluid are necessary to achieve this outcome. 
Hence the primordial state of matter is described as a 
strongly-coupled quark gluon plasma (sQGP). 

Quantum chromodynamics (QCD) produces the bulk 
of the mass of normal matter. At T = it is characterised 
by confinement and dynamical chiral symmetry breaking 
(DCSB), phenomena that are represented by a range of 
order parameters which all vanish in the sQGP. Under- 
standing the sQGP therefore requires elucidation of the 
behaviour and properties of quarks and gluons within 
this state. Perturbative techniques have been developed 
for use far above T c ; viz., the hard thermal loop (HTL) 
expansion [4, 5], which has enabled the computation of 
gluon and quark thermal masses tut ~ gT and damping 
rates 7t ~ g 2 T, with g = g(T) being the strong running 
coupling. It also suggests the existence of a collective 
plasmino or "abnormal" branch to the dressed-quark dis- 
persion relation, which is characterised by antiparticle- 
like evolution at small momenta [6]. 

Owing to asymptotic freedom, the running coupling in 
QCD increases as T — >• T c + . Therefore, a simple interpre- 
tation of the HTL results suggests the plasmino should 
disappear before T c is reached because increases more 
rapidly than tut and jr/mr ~ 1 invalidates a quasi- 



particle picture. On the other hand, lattice-regularised 
quenched- QCD suggests that the plasmino branch per- 
sists in the vicinity of T c [7]. It is necessary to resolve 
the active degrees of freedom in the neighbourhood of 
T c because the spectral properties of the dressed-quark 
propagator are intimately linked with light-quark con- 
finement [8] and it is the long-range modes which might 
produce strong correlations. 

When addressing issues concerning the dressed-quark 
propagator it is natural to employ the gap equa- 
tion, which is one of QCD's Dyson- Schwinger equations 
(DSEs) [9-11]. Equations of this type are ubiquitous 
in physics and in QCD the DSEs are distinguished by 
their ability to unify the analysis of confinement and 
DCSB within a single nonperturbative, Poincare covari- 
ant framework. Extensive work within this approach has 
shown that below T c dressed-gluons and -quarks are con- 
fined and chiral symmetry is dynamically broken in the 
chiral limit [10, 11], and that deconfinement and chiral 
symmetry restoration occur via coincident second-order 
phase transitions [12-14]. Herein we use the DSEs to 
elucidate the active fermion quasiparticles for T >T C . 

On the domain T > T c the chiral-limit dressed-quark 
propagator can be written 

S(iu n ,p) -ij • pa A (^n,p 2 ) ~ ij4^nCrc(^n,P 2 ) , (1) 

where uj n = (2n + 1)ttT, n G Z, is the fermion Matsubara 
frequency. There is no Dirac-scalar part because chiral 
symmetry is realised in the Wigner mode. The retarded 
real-time propagator is found by analytic continuation 

S R (w,p) = S{iu n ,p)\ iuJn ^ +iv+ (2) 

and from this one obtains the spectral density 

p(u J ,p) = -2SS R (u J ,p). (3) 

Equations (2) and (3) are equivalent to the statement: 
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Notably, if one requires a nonnegative spectral density, 
then Eq. (4) is only valid for T > T c ; i.e., on the decon- 
fined domain [12]. 

For an unconfined dressed-quark propagator of the 
form in Eq. (1), the spectral density can be expressed 



-m n 



p(u,p) = p+{uj,p 2 )P+ + p-(uj,p 2 )P- , 



(5) 



where P± = (74 ± • u p )/2, u p • p = \p\, are operators 
which project onto spinors with a positive or negative 
value for the ratio T~L := helicity/chirality: T-L = 1 for a 
free positive-energy fermion. The spectral density is in- 
teresting and expressive because it reveals the manner by 
which interactions distribute the single-particle spectral 
strength over momentum modes; and the behaviour at 
T/0 shows how that is altered by a heat bath. As with 
many useful quantities, however, it is nontrivial to eval- 
uate p(lj, \p\). Nonetheless, if one has at hand a precise 
numerical determination of the dressed-quark propaga- 
tor in Eq. (1), then it is possible to obtain the spectral 
density via the maximum entropy method (MEM) [15]. 

We obtain the chiral-limit dressed-quark propagator 
from the gap equation 

S(iuj n ,p)~ l = Zztf • p + Z 2 ij4^n + H'(iu n ,p) , (6) 
H'(iuj n ,p) = ij ■pT,' A {iuj n ,p) + ij4UJ n E' c (iu n ,P) 

/d?q 
J2^9 2 D^ (u n -u u p-q) 

x^tT B V F JfiS(iuJi,q)T u (uj n ,uJi,p,q), (8) 

where: Va = -Zftf-p/p 2 , Vc = -Z\i^/u n \ D^ v is the 
dressed-gluon propagator; T v is the dressed-quark-gluon 
vertex; and ^1,2, Z^ 2 are > respectively, the vertex and 
quark wave function renormalisation constants [12, 16]. 

The gap equation is determined once the kernel is spec- 
ified. Herein we work at leading-order in the symmetry- 
preserving truncation scheme of Ref. [17] and employ a 
phenomenologically-efflcacious one- loop renormalization- 
group-improved interaction [16]. Namely: 

g 2 D fIU (u n - LOup- q)T v (u n ,ui,p,q) 
= [P^(kn)D T (k n ) + P£ u (kn)D L (k n )h^ (9) 



where kn := (ft, k) = (cj n - uj^p- q)\ 



0, 



P^(k n ) 



p and/or v = 4 , 



kj k q 
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with P» v + P£ u = 5^ - k&k&fk^ and 



D T {k n )=V{kl,0), D L (kn)=-D(k 2 Q , ml), (11) 
V(klm 2 g )=4n 2 D S ^e-^ 



ln[r+(l + ^/A 2 QCD )2] 



(12) 



with J'(sn) = (l-exp(-sn/4m|)/sn, sq = ft 2 +^ 

r = e 2 - 1, ra t = 0.5 GeV, j m = 12/25, and Aq^ 4 = 
0.234 GeV. For pseudoscalar and vector mesons with 
masses < 1 GeV, this interaction provides a uniformly 
good description of their T = properties [18] when 
u = 0.4 GeV, D = (0.96 GeV) 2 . In generalizing to T ^ 0, 
we have followed perturbation theory and included a De- 
bye mass in the longitudinal part of the gluon propagator: 
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(16/5)T 2 . 



A justification of the kernel is readily provided. At 
T = it reproduces the results of perturbative QCD for 
p 2 > 2 GeV 2 , so any model-dependence appears only in 
the infrared, and provides a unified description of light- 
vector and -pseudoscalar mesons. It also predicts a mo- 
mentum dependence for the dressed-quark propagator 
that is qualitatively in agreement with results from nu- 
merical simulations of lattice-QCD [19]. The extension to 
T > preserves the agreement with perturbative QCD 
at large spacelike momenta. Finally, in employing the 
kernel we obtain coincident second-order deconfinement 
and chiral symmetry restoring transitions for two mass- 
less flavors at T c = 0.14 GeV, which is 10% smaller than 
that obtained in Ref. [20]. 

One insufficiency of the interaction defined above is 
that D, the parameter expressing its infrared strength, is 
assumed to be T-independent. Since the nonperturbative 
part of the interaction should be screened for T > T c , we 
remedy that by writing D — >> D(T) with 



D(T) 



D, 



T<T D 



MT/Kqcd] ' T ~ Tp 



(13) 



where T p is a "persistence" temperature; i.e., a scale be- 
low which nonperturbative effects associated with con- 
finement and dynamical chiral symmetry breaking are 
not materially influenced by thermal screening. Loga- 
rithmic screening is typical of QCD and with a = 0.028, 
b = 0.56 our numerical solutions yield vcit = 0.8 T for 
T > 2T C ; viz., a thermal quark mass consistent with 
lattice-QCD [7]. We usually take T p = T c herein. 

We have computed the spectral density by employing 
the MEM in connection with the solution of our gap equa- 
tion. Notably, the behaviour changes qualitatively at T c . 
Indeed, employing a straightforward generalisation of the 
inflexion point criterion introduced in Refs. [11, 21], one 
can readily determine that reflection positivity is violated 
for T <T C . This signals confinement. On the other hand, 
the spectral function is nonnegative for T > T c . 

In Fig. 1 we depict the T > T c -dependence of the loca- 
tions of the poles in p(u),p = 0); i.e., the thermal masses. 
We anticipated that spectral strength would be located at 
uj+(p = 0) and 0J-(p = 0) = — oj+(p = 0), corresponding 
to the fermion's normal and plasmino modes at nonzero 
temperature. However, it is striking that on a measur- 
able T-domain, spectral strength is also associated with 
a quasiparticle excitation described by uoq(p = 0) = 0. 
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FIG. 1. (Color online) Left panel - Temperature-dependence 
of the dressed-quark thermal masses. Notably, spectral 
strength is associated with a massless-mode. Right panel - 
T-dependence of the residue associated with that zero mode. 

The appearance of this zero mode is an essentially non- 
perturbative effect. It is an outgrowth of the evolution 
in-medium of the gap equation's T = Wigner-mode 
solution and analogous to this solution's persistence at 
nonzero current-quark mass in vacuum [22]. 

The spectral density possesses support associated with 
this zero mode on T G [0,T S ]. In fact: all the Wigner- 
phase spectral strength is located within this mode at 
T = 0; it is the dominant contribution for T > T c ; and, 
while it is dominant, it is the system's longest wavelength 
collective mode. On the other hand, as evident in the 
right panel of Fig. 1, the mode's spectral strength dimin- 
ishes uniformly with increasing T and finally vanishes at 
T s « 1.35 T c . Then, for T > T s the quark's normal and 
plasmino modes exhibit behavior that is broadly consis- 
tent with HTL calculations. This is apparent in Fig. 1 
and in a comparison between the upper and lower panels 
of Fig. 2. Given these observations, we judge that the 
system should be considered a sQGP for T G [T C ,T S ], 
whereupon it contains a long-range collective mode. 

We observe that the HTL approach is perturbative and 
only applicable for T ^> T c . Hence it could not have 
predicted the zero mode's existence. Numerical simula- 
tions of lattice-QCD, on the other hand, are nonpertur- 
bative. However, it is practically impossible in contempo- 
rary computations to exactly preserve chiral symmetry. 
This can plausibly explain the absence of the zero mode 
in lattice simulations because any source of explicit chiral 
symmetry breaking heavily suppresses the mode [22]. 

The upper-left panel of Fig. 2 depicts the dispersion 
relations for all dressed-quark modes that exist for T < 
T s and the behavior of their associated residues. On this 
sQGP domain the dispersion relations are atypical, with 

/ M xp~o -0.2 \p\, , . 

^(IpI) = ^ +o.3 1^|, ( 14 ) 

u;o(|p1) P = 0.80|p1. (15) 

Notwithstanding this, all realise free-particle behaviour 
for \p\ ^> T. We note and emphasise that the usual 




FIG. 2. (Color online) Upper panel, left - Quasiparticle dis- 
persion relations, cj±,o(p), at temperature T = 1.1T C . The 
diagonal dotted line is the free-fermion dispersion relation at 
this T. Upper panel, right - The residues associated with 
these quasiparticle poles. Lower panel - Same information 
for T = 3T C , whereat the zero mode has vanished. 

spectral sum rules are satisfied. Indeed, the identity 

^ — u ' p± (u,',\P\) = Zf\0l (16) 

assists in understanding the momentum-dependence in 
the upper- left panel of Fig. 2. The upper-right panel dis- 
plays the momentum-dependence of the pole residues: 
spectral support is located completely in the normal 
mode for \p\ ^> T; i.e., on the perturbative domain. 

The lower panel of Fig. 2 characterises the behaviour 
of p(cj, |p|) at on the high-T domain. In agreement with 
HTL analysis [5], expected to be valid thereupon, we find 
only normal and plasmino modes, with 

cu±(\p\) P =m T ±0.33\p\. (17) 

The plasmino dispersion law exhibits the expected min- 
imum, in this case at \p\/T ~ 1; and both uj±(\p\) ap- 
proach free-particle behaviour at |p| ^> T, with that of 
the plasmino approaching this limit from below. The 
lower-right panel shows that the contribution to the spec- 
tral density from the plasmino is strongly damped and 
contributes little for p > 2T. These results are in-line 
with those obtained via simulations of lattice-QCD [7]. 

Equation (13) is a model and it is natural to enquire 
after its influence. None of our results are qualitatively 
altered by varying T p but, as one would expect, the width 



of the sQGP domain expands slowly with increasing T p ; 
e.g., a 50% increase in T p produces a 30% increase in T s . 

Whilst we used the MEM to compute p(u), \p\) from the 
completely self-consistent numerical solution of a rainbow 
truncation of QCD's gap equation, the appearance of a 
third and long-wavelength mode in the dressed-fermion 
spectral density on a material temperature domain above 
T c has also been observed in one-loop computations of the 
fermion self-energy, irrespective of the nature of the bo- 
son which dresses the fermion [23]. This mode appears 
for T > mo-, where mo is an infrared mass-scale associ- 
ated with the boson. In our case, mo = 0.12 GeV. Where 
a comparison is possible, the dependence of our spectral 
density on (lj, \p\,T) is similar to that seen in the one- 
loop analyses of model gap equations. In analogy with a 
similar effect in high-temperature superconductivity [25] , 
that behaviour has been attributed to Landau damping, 
an interference phenomenon known from plasma physics. 

Notably, we find that a coupling to meson-like correla- 
tions in the gap equation is not a precondition for the ap- 
pearance of the zero mode because such correlations are 
absent in the rainbow truncation [13]. On the other hand 
our gap equation's kernel is characterised by an interac- 
tion that features an infrared mass-scale mo < T c and 
supports dynamical chiral symmetry breaking at T = 0. 
We anticipate that the zero mode will markedly affect 
colour-singlet vacuum polarisations on T G [T C ,T S ]. This 
could be explicated using the methods of Refs. [26]. 

Hadron physics experiment has presented us with the 
fascinating possibility that a near-perfect fluid might 
have been a key platform in the universe's evolution. We 
have sought to provide new insights into this possibility 
by employing the Dyson- Schwinger equations (DSEs), a 
tool used widely in many branches of physics and which 
in QCD provides a nonperturbative, continuum tool that 
unifies the treatment of confinement, dynamical chiral 
symmetry breaking and observable phenomena. 

In pursuing this aim, we solved self-consistently a rain- 
bow truncation of the gap equation for massless 2-flavour 
QCD, employing a kernel whose temperature-dependence 
is constrained by T = hadron physics phenomenology 
and T / lattice-QCD results, and found that chiral- 
symmetry restoration and deconfinement occur together 
at a temperature T c = 140 MeV. We subsequently com- 
puted the quark spectral density from the gap equa- 
tion's solution for T > T c using the maximum entropy 
method, thereby demonstrating its potential when based 
on accurate input. Remarkably, on a significant domain 
T/T c e [1,1 + A], A - 0.4 - 0.8, the self-consistently 
determined Wigner phase supports a zero mode, despite 
the absence of meson-like correlations in our gap equa- 
tion's kernel. This mode contains the bulk of the spectral 
strength for T > T c and so long as this mode persists, 
the system may reasonably be described as a strongly- 
interacting state of matter. If, as we argued, the existence 



of this long- wavelength mode is model-independent, then 
it is natural to anticipate that a strongly-interacting state 
of matter should precede the QCD phase transition. 
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